Stable Sarma State in Two-band Fermi Systems 
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We investigate fermionic superconductivity with mismatched Fermi surfaces in a general two- 
band system. The exchange interaction between the two bands changes significantly the stability 
structure of the pairing states. The Sarma state with two gapless Fermi surfaces which is always 
unstable in single-band systems, can be the stable ground state in two-band systems. To realize a 
visible mismatch window for the stable Sarma state, two conditions should be satisfied: a nonzero 
inter-band exchange interaction and a large asymmetry between the two bands. 
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I. INTRODUCTION 



The Cooper pairing with mismatched Fermi surfaces, 
which has been investigated many years ago^"^, promoted 
new interest in the study of new superconducting materi- 
als in strong magnetic field and ultracold fermions due to 
the realization of superfluidity in resonantly interacting 
Fermi gases. The well-known theoretical result for s-wave 
weak coupling superconductors is that, at a critical mis- 
match, called Chandrasekhar-Clogston limit (CC limit) 
he = 0.707Ao where Aq is the zero temperature gap, a 
first order phase transition from the gapped BCS state 
to the normal state occurs^. Further theoretical stud- 
ies showed that the inhomogeneous Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state^ may survive in a narrow win- 
dow between he and /ifflo = 0.754Ao. However, since 
the thermodynamic critical field is much smaller than 
the CC limit due to strong orbit effecli^, it is hard to 
observe the CC limit and the FFLO state in ordinary 
superconductors. In recent years, some experimental ev- 
idence for the FFLO state in heavy fermion supercon- 
ductor^, high temperature supercondutors^ and organic 
superconductors^ were found. More recently, in the study 
of ultracold atoms, Fermi superfluidity with population 
imbalance was realized by MIT and Rice groups indc- 
pendently^. The ultracold fermion experiments has pro- 
moted a lot of theoretical work o^i^i^°'^^ on the superfluid- 
ity mechanism and the phase diagrams for the crossover 
from Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein 
Condensation(BEC) The problem of imbalanced pair- 
ing is also related to the study of color superconductivity 
and pion superfluidity in dense quark matter— i^^. 

While most of the theoretical works focus on the in- 
homogeneous FFLO state, we in this paper are inter- 
ested in the homogeneous and gapless Sarma state^. For 
weak coupling superconductors, the Sarma state is lo- 
cated at the maximum of the thermodynamic potential 
of the system, and therefore can not be the stable ground 
state. This was called Sarma instability many years agoA. 
The thermodynamic instability of the Sarma state can 
be traced to the existence of gapless fermion excitations 
which cause a very large density of state at the gap- 
less Fermi surfacesi*^. To realize a stable Sarma state, 
one should have some mechanism to cure the instabil- 



ity. Forbes et.al>i^ proposed that, a stable Sarma state 
is possible in a model with finite range interaction where 
the momentum dependence of the pairing gap cures the 
instability. On the other hand, when the attractive inter- 
action becomes strong enough which can be realized in 
ultracold fermion experiments, the stability of the Sarma 
state can be changed. While the homogeneous Sarma 
state is always unstable at the BCS side of the BCS- 
BEC crossover, it becomes stable in the deep BEC re- 
gion^. However, this stable Sarma state at the BEC side 
is not the original "interior gap" or "breached pairing" 
state with two gapless Fermi surfaces proposed by Liu 
and Wilczek^^. Since the fermion chemical potential be- 
comes negative in the BEC region, the Sarma state in 
this case possesses only one gapless Fermi surface, and 
the matter behaves like a Bose-Fermi mixtureii. 

In this paper, we focus on how the multi-band struc- 
ture which may be realized in solid materials and optical 
lattices changes the stability of the Sarma state. We con- 
sider a general two-band Fermi system, and show that the 
inter-band exchange interaction can cure the Sarma in- 
stability and the Sarma state can be the stable ground 
state in visible parameter regions. 

The multi-band theory of BCS superconductivity was 
firstly introduced by Suhl et a\M in 1959 to describe the 
possible multiple band crossings at the Fermi surface. 
The two-band model has been applied to the study of 
high- Tc superconductors^^ to effectively describe the par- 
ticular crystalline and electronic structure. Recently, it 
is found that, the material MgB2 is a standard two-band 
superconductor— and many experimental data can be 
explained by the two-band model of BCS superconduc- 
tivity. Multi-band Fermi systems may be realized experi- 
mentally with ultracold atoms in optical lattice^. For ex- 
ample, if we confine the cold atoms in a one dimensional 
periodic external potential, the band structure will form 
in the confined direction, and the matter can be regarded 
as a multi-band system in two-dimensions. In this case, 
by adjusting the coupling strength, one can study the 
possible BCS-BEC crossover in multi-band systemai^i^^. 
The inter-band physics in optical lattices is recently stud- 
ied24, and the multi-gap superfluidity is also possible in 
nuclear matter—. 

The paper is organized as follows. In Section [H] we 
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give an introduction to the Sarma state in single-band 
systems. We discuss the stabihty of the Sarma state in a 
general two-band model in Section IIIII and summarize in 
Section HVl 

II. SARMA STATE IN SINGLE-BAND MODEL 



Whether the Zeeman energy imbalance h or the spin den- 
sity imbalance S is experimentally adjusted depends on 
detailed systems. For cold atoms, the spin density im- 
balance S is directly tuned, but in superconductors, the 
Zeeman splitting h is adjusted via an external magnetic 
field. 



Before discussing the stability of Sarma state in two- 
band systems, we in this section give a brief introduction 
to the Sarma state in single-band systems. We start from 
the following standard BCS-type Hamiltonian 

2m 



H 



(1) 



We constrain ourselves to discuss systems at zero temper- 
ature where the BCS mean field theory can be applied 
even at strong coupling—. In the mean field approxima- 
tion, the Hamiltonian is approximated by 



mi 



+ $(r)V|(r)Vl(r)-HH.c. 



u 



(2) 



where $(r) = —U{ip^{r)ijj^{r)) is the order parameter 
field of superconductivity. For homogeneous supercon- 
ductivity, the thermodynamic potential f2 can be ob- 
tained by using the standard diagonal method-'^®. It can 
be expressed as 



n 



A2 

U 



(3) 



with the definition of energy dispersions = k;^/(2m) 



fi, El, = y^fTA^", eI = Ei, + h and E^ = E^-h, where 
H = (/i-i -I- /X|)/2 and h = (/if — ^J-l)/'2, are, respectively, 
the averaged and mismatched chemical potentials, and 
A is the modulus of ^(r). 

Without loss of generality, we set h > 0. The possible 
ground state of the system corresponds to the stationary 
point of the thermodynamic potential fl. This gives the 
so-called gap equation 




A = 0. 



(4) 



To properly achieve strong coupling, the chemical poten- 
tials should be renormalized by the number equations. 
The number density n and spin density imbalance S can 
be evaluated as 



(27r)3 



1 



El, 



(5) 



A. Stability Analysis 

If a solution of the gap equation is the ground state of 
the system, it should be the global minimum of the ther- 
modynamic potential fl^^-^^. The condition for a local 
minimum of f2 is that 



dA 



= 0, 



9A2 



> 0. 



(6) 



The first condition corresponds to the gap equation and 
the second order derivative / — d^ft{A) / dA^ can be eval- 
uated as 



1 = 



d^k A2 \e{El) 



(2^)3 Ei 



(7) 



Let us study the Sarma state with h > A which induces 
a nonzero spin density imbalance 5. At weak coupling, / 
can be approximately evaluated as 



•\/2m/i 



1 



h<d{h 



Vh^ - A2 



(8) 



which shows that d'^^l{A) /dA^ is always negative and 
therefore the Sarma state is unstable. 

To achieve the BCS-BEC crossover, we renormalize the 
coupling constant with the two-body scattering length a^. 



47ra, 



1 

U 



d^k 



(27r)3 k2 ■ 



(9) 



In this case, we first solve the coupled gap and number 
equations at fixed total density n = fcp/(37r2). The re- 
sult can be expressed^ as a function of the dimensionless 
coupling parameter g = 1/ (kpag) and the population im- 
balance P = 5/n. The numerical calculations show that, 
the key quantity / is always negative at the BCS side 
of the resonance(a5 < 0, g < 0) where the Sarma state 
has two gapless Fermi surfaces, but the Sarma state can 
be a stable ground state in the strong coupHng BEC re- 
gion (roughly for g > 2.2) where the chemical potential 
fi become negative. However, this Sarma state has only 
one gapless Fermi surface and is different from the Fermi 
surface topology of the so-called breached pairing state. 



B. Solution at Weak Coupling 

At weak coupling where the chemical potential fi is well 
approximated by the Fermi energy ep, the gap equation 
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can be approximated by 



III. SARMA STATE IN TWO-BAND MODEL 



UN- I 



A = 0, 



(10) 



where N is the density of state for each spin state at the 
Fermi surface, and A is the energy cutoff which plays the 
role of Debye energy hujjj in solids. After the integration 
and using the condition A ^ A, wc find 



1 2A 
— -ln- + e(/.-A)ln 



h + Vh^- A2 
A 



A = 0. 



(11) 

There are three possible solutions to the gap equation 
(fTTjl for h ^ Q. The first is the trivial normal phase with 
An — 0. The second corresponds to the ordinary fully 
gapped BCS solution satisfying A > ft,, 



ABcs = Ao = 2Ae-i/(^^). 



(12) 



The third solution, i.e., the gapless Sarma state satisfying 
A < h, can be analytically evaluated via the comparing 
with the BCS solution. It is2^ 



As = VAo(2/i- Ao). 



(13) 



Using the weak coupling approximation, the grand po- 
tential Q for various solutions can be expressed as 



A2 

n = — 

u 



2N 



+ i^/e + A2 - h)eih - + A2)j . (14) 

Performing the integral over ^, and using the condition 
A <C A as well as the gap equation to cancel the cutoff 
dependence, we have 



n 



N 



e{h- A)Nhy/h? 



(15) 



Note that we have set the grand potential of the normal 
state at ft = to be zero, il^{h = 0) = 0. To see why the 
Sarma state is always thermodynamically unstable, one 
should calculate the grand potential differences between 
Sarma and other two states^^, 



ns-n^ = y(Ao-2ft)2, 



(16) 



which confirm that the Sarma state always has higher 
potential than the BCS and normal states. As a conse- 
quence, there exists a first order phase transition from 
BCS state to normal state. From the result 



BCS 



N 



(2ft; 



(17) 



the transition occurs at the CC limit of BCS supercon- 
ductivity, he = Ao/\/2- 



We in this section turn to the two-band model. Since 
the goal of this paper is to search for the possibility of sta- 
ble Sarma state in general two-band Fermi systems, we 
consider a continuum Hamiltonian and neglect the details 
of the band structure in different systems. We will show 
that the key point is the inter-band scattering which can 
make the Sarma state stable in multi-band systems. The 
possible complicated lattice structure in various materi- 
als and optical lattices will not qualitatively change our 
conclusion. The obtained conclusion is generic and may 
be useful for the study of superconducting materials and 
ultracold atom gases. 

The continuum Hamiltonian of the two-band model 
can be written as^^ 



H 



2m„ 



(18) 



where I'jX = 1,2 denote the band and a =t,i the di- 
rection of fermion spin. In superconductors, the band 
degrees of freedom usually come from the particular crys- 
talline and electronic structure of the materials. In ul- 
tracold atom gases, these degrees of freedom may come 
from different hyperfine states or different atom species 
or the external periodic lattice potential. In general case, 
the effective fermion mass depends only on the band in- 
dex, but the chemical potential is related to both the 
band and spin indexes due to the existence of external 
magnetic field or population imbalance. The constants 
Uii = Ui and U22 = U2 are the intra-band couplings, 
and U12 = U21 = J is the inter-band exchange coupling. 
For vanishing J, the model is reduced to a simple system 
with two independent bands. In the following we focus 
on how the inter-band coupling J changes the stability 
of the Sarma state. 

We first calculate the thermodynamic potential of the 
two-band Hamiltonian. In the mean field approximation, 
the Hamiltonian is approximated by 



+ ^[c/2|$i(r)|2 + t/i|$2(r)P 
- J($t(r)<i>2(r)-K<i>;(r)ci>i(r)) }, 



(19) 



where $i/(r) — — X^a ^iyA(^Ai('")V'AT('")) ^^'^ ^^'^ order 
parameter fields of the superconductivity, and G is de- 
fined as G = U1U2 — ■ For homogeneous superconduc- 
tivity, the thermodynamic potential il of this two-band 
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model can be obtained by using the standard diagonal 
methodiS. It can be expressed as 



n ^ - [U2AI + UiAj - 2 JA1A2 cos{ipi - (p2)] (20) 



E 



(2^ 



with the definition of energy dispersions ^ki^ = 
kV(2m,) - fi,, E^, = V?^, + A2, = E^, + /i, 

and e]^^ = E'ki/ - K, where = (/Xi,| + [ivi)!"! and 
^1/ = (Mi^t ^ are, respectively, the averaged and 

mismatched chemical potentials, and Ay the modulus 
of and iy9y their phases through the definition '^y — 
AyC^^" . Without loss of generality, we take hy > 0. For 
J > 0, the choice of (pi = (p2 is favored, otherwise there 
is (fii — Lp2 + TT. We assume J > and set ipi ~ ip2- 

The possible ground state of the system corresponds 
to the stationary point of the thermodynamic potential 
f2. This gives the so-called gap equations 



G 



(2^)3 2£;k. 



Ay 



Ay^O 



(21) 



with 9—1 for v — 2 and 9 = 2 for v — 1. The gap 
equations are essentially the same as derived ini^. To 
properly achieve strong coupling, the chemical potentials 
should be renormalized by the number equations. The 
number equations for the fermion density Uy and density 
imbalance Sy for the v-th band can be evaluated as 



Sy = riy-l — = 



(2^ 
(2^ 



1 - p^mL) 



^i~Ety), 



(22) 



and the total density n and total density imbalance S of 
the system are defined as n = ni + 712 and S = Si + 62- 



A. Stability Analysis 

Let us now discuss qualitatively what happens when 
the mismatch hy increases. For vanishing mismatch, the 
system is in a fully gapped BCS state with Ai = Aio 
and A2 = A20, and the spin density imbalance S is zero. 
With increasing hy, while the BCS state is still a solution 
of the gap equations, there may appear another solution 
(Sarma) where at least one of the pairing gap Ay is less 
than the corresponding mismatch, namely hy > Ay. In 
this state, the dispersion of the quasi-particle E^^^ be- 
comes gapless and the system has a nonzero spin density 
imbalance S. Note that the normal state with vanishing 
condensate is always a solution of the gap equations and 
becomes the ground state when both hi and are large 
enough. 

Different from the conventional Sarma state in single- 
band models, we may have two types of Sarma states in 



two-band systems. The first type (type I) is the solution 
where both mismatches are larger than the correspond- 
ing pairing gaps, namely hi > Ai and /i2 > A2. For 
this type, there exist gapless excitations in both bands. 
The second type (type II) is the solution where only one 
mismatch is larger than the corresponding pairing gap, 
hi > Ai and /12 < A2 or hi < Ai and /12 > A2. For 
this type, gapless excitations exist only in one band. We 
will show in the following that the stabilities of these two 
types of Sarma states are quite different. 

A numerical example which supports the above argu- 
ment is shown in FigU] for two symmetric bands with 
Ui — U2- For the sake of simplicity, in our numer- 
ical examples presented here, we assume the same ef- 
fective masses, chemical potentials and mismatches for 
the two bands, i.e., mi = TO2 = m, ^1 = /i2 = /i 
and hi = h2 = h, this means that only the total den- 
sity n and total spin density imbalance S can be ad- 
justed. We write Ui and U2 in terms of the s-wave scat- 
tering length Qy with a momentum cutoff fco, Uy^ = 
— m/(47ray) + J^^^^^^,^d^'k/{2TT)^m/'k'^. Our qualitative 
conclusions do not depend on the used regularization 
scheme. In the case of Ui — U2, the solutions of the gap 
equations are distributed symmetrically in the Ai — A2 
plane. Besides the familiar BCS and normal states which 
are, respectively, the global minimum and a local mini- 
mum in Figlll we have some Sarma states in the potential 
contour. The Sarma states C and D are of type I, and 
C is the global maximum and D indicates two saddle 
points. The type II Sarma states are marked by A and 
B, corresponding, respectively, to two local minima and 
two saddle points. 

If a solution of the gap equations is the ground state of 
the system, it should be the global minimum of the ther- 
modynamic potential fii^i^. The condition for a local 
minimum of is that the matrix 



M 



3^o(Ai A2) a^n(Ai,A2) 

aA^ OAiSAs 

a^tXAi.Aa) 3^0(Ai,A2) 
aAaSAi dM 



(23) 



should have only positive eigenvalues, namely detM > 
and TrAi > 0. The second order derivatives can be 
evaluated as 



a2rj(Ai,A2) _ 2jAy 



dAl 



G Ay 



ly 



92f](Ai,A2) _ 2J 
dAydAy ^ ^~G 

with the quantities ly defined as 



d^k Al 



Eky 



SiEty) 



(24) 



(25) 



For vanishing inter-band coupling ,7 = 0, the stability 
condition becomes 



/l/2>0, /i+/2>0. 



(26) 
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FIG. 1: The thermodynamic potential contour f2(Ai,A2) 
for two symmetric bands with Ui = U2- A proper unit 
is chosen such that the Fermi energy ep — 200. The val- 
ues of the other parameters are 

with Uo = 47r/(mfcF), (fcpai)"^ = {kFa2)~'^ — —0.5 and 
h — 75, where kp ~ ^/2mfj, is the Fermi momentum. The 
band on the right shows the relative strength of Q corre- 
sponding to different colors. For the parameter setting, we 
have f/i = (72 ~ 0.0156f/o, and hence J ^Ui,U2. 



Note that the properties of the functions /i and I2 are 
the same as the function / defined in the last section. 
Thus at the BCS side, namely for ai < and a2 < 0, the 
Sarma state is unstable. 

Now we discuss how the inter-band coupling J modifies 
the Sarma instability at the BCS side. For J 7^ 0, the 
stability condition reads 

For the type I Sarma state with hi > Ai and /12 > 
A2, we have Ii < and I2 < 0. In this case, we can 
exactly prove that the above two inequalities can not 
be satisfied simultaneously. This type of Sarma state 
should correspond to the maximum or saddle point of the 
thermodynamic potential and is hence unstable, like the 
points C and D in Fig[TJ However, the situation changes 
for the type II Sarma state. Without loss of generality, 
let us discuss the case with hi > Ai and ^2 < A2. In 
this case, only the first band is gapless, and hence Ii < 
and I2 > 0- From I2 > 0, the above two inequalities are 
equivalent to the following one 

Even though Ii < 0, this condition can be satisfied, pro- 
vided that a nonzero inter-band coupling J is turned 



on. Suppose the solution of the gap equations satisfies 
Ai ^ A2 and Ai is not quite close to hi, which cor- 
responds to the case with large polarization 61, the first 
term in (|28p is large but the modulus of the second term 
is relatively small, and therefore the stability condition 
can be satisfied, like the point A at the up-left cornel in 
FiglH However, on the other hand, for Ai < hi which 
corresponds to the case with small polarization 61 — s- 0, 
the absolute value of /i is very large, and the Sarma state 
maybe unstable, which corresponds to the saddle point 
B in the upper part of Fig[TJ 

We conclude that, in two-band Fermi systems with 
non-zero inter-band pairing interaction, the Sarma state 
can become at least the local minimum of the thermo- 
dynamic potential, and therefore should be a potential 
candidate of the ground state. 

However, the condition J 7^ is not sufficient for us 
to have a real stable Sarma state. For the case with 
two symmetric bands shown in FiglU we found that the 
global minimum is always the BCS or normal state for 
any mismatch h, which means that the Sarma state can 
not be the ground state even though it can be a local 
minimum. However, this can be significantly changed if 
some asymmetry between the two bands, such as unequal 
couphngs Ui ^ U2, is turned on. In Figl2]and Fig[31 we 
show the potential contour with Ui ^ U2 for three val- 
ues of h. In this case, the number of Sarma solutions is 
largely suppressed due to the asymmetry, especially the 
state C in Fig[l]as the global maximum of disappears. 
Without regard to the saddle points which are impossi- 
ble to be stable solutions, the only Sarma state marked 
in Figl2] and Fig[3] appears to be the global minimum of 
the system, when the mismatch /i is in a suitable region. 
From the top to the bottom in Figl^] and FiglJl when 
the mismatch h increases, the global minimum changes 
from the BCS state to the Sarma state and then to the 
normal state. In contrast to the conventional single band 
model where only one first order phase transition from 
the BCS to normal state is predicted, we have in this 
two band system two first order phase transitions when 
h increases: The first is from the BCS to Sarma state, 
and the second is from the Sarma to normal state. The 
first order phase transition from BCS to Sarma state was 
found in^^ by considering the momentum structure of 
the pairing gap. For ultracold atom gases, the chemi- 
cal potential mismatch h should be replaced by the spin 
population imbalance S. However, the phase structure 
should be essentially independent of the assembles one 
usedi^i^, we here do not consider the case with fixed S. 

Let us compare the numerical results presented in Fig [2] 
and FigO In Figl3l the coupling asymmetry is much 
larger than that in Figl^l we have A10/A20 — 1.5 in FigO 
and Aio /A20 ~ 4 in FiglJl We find that the h window for 
the Sarma state is wider when the coupling asymmetry 
becomes larger. In Figl2l the window for Sarma state is 
roughly from /i = 52 to ft, = 71, and the CC limit is about 
he ~ A20. In Fig|31 this window is roughly from h — 20 
to ft = 70, and the CC limit is about he — 2.7A2o. 



20 40 60 80 100 120 



20 40 60 30 100 120 



FIG. 2: The thermodynamic potential contour n(Ai, A2) for 
two different bands with (ftpffli)"^ = —0.5 and (fcpa2) 
-0.8 at /i = 45 (top), 60 (middle) and 75 (bottom), 
other parameters are the same as that in FigIT] 



FIG. 3: The thermodynamic potential contour r2(Ai, A2) for 
two different bands with (fcpai)"^ = —0.5 and (A;Fa2)~^ = 
-f.5 at ft = 20 (top), 25 (middle) and 70 (bottom). The 
other parameters are the same as that in Fig[T] 



B. Solutions at Weak Coupling 



and express the gap equations of our two band model in 
terms of it, 



At weak coupling, the same tricks used in Section |TT] 
can be employed. For convenience, we define here a func- 
tion 



F(A,/i) = 







(29) 



ln^-e(/i- A)ln- 



V/i2 - A2 



Al —A2 = 0, 

GNi 



A, 



J 

cm 



■Al = 0, (30) 



where Ni and N2 are the densities of state at the Fermi 
surfaces for the two bands. Unlike the single band model, 
the above coupled gap equations can not be solved ana- 
lytically. With the numerical solutions Ai and A2, the 
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thermodynamic potential can be evaluated as 



2 " 



.(31) 



For the sake of simplicity, we consider the case hi = 
h2 = h which corresponds to the realistic two-band su- 
perconductors in strong magnetic field. Let us assume 
UiNi > U2N2 which leads to Ai > A2. According 
to the stability analysis, we have three possible ground 
states: l)The normal state with Ai = A2 = 0; 2)Thc 
gapped BCS state with energy gaps Ai = Aio > h and 
A2 = A20 > h; 3)The gapless Sarma state where only 
Ai > h but A2 < h. We focus here on the case with 
A2 <C Ai and J <C \/UiU2- In this case, the solution of 
Ai is approximately independent of h and is given by 



Ai 



^10 



(32) 



and the Sarma solution for A2 is determined by the fol- 
lowing equation 



In- 



A. 



20 



A2 
^2 



JAio 

U1U2N2 



1 



where A20 is obtained by the equation 



U2N2 



In 



2A 
A~ 



J A 



1 

A^ 



10 



20 



U1U2N2 A20' 



(33) 



(34) 



We have numerically checked that the above approxi- 
mation is sufficiently good for the scaled solution y = 
A2/A20 as a function of a; = /1/A20. Note that for J = 
the conventional Sarma solution y = \j2x — 1(0.5 < a; < 
1) is recovered, but for J ^ the Sarma solution is qual- 
itatively changed: y — ^ can not be a solution and there 
exist solutions for a; > 1. The solutions for both cases of 
J = and J ^ are illustrated in FiglH We find that 
for J ^ the Sarma solution is quite different from the 
conventional result. Unlike the well-known Sarma solu- 
tion which exists in the region 0.5 < a; < 1, for J ^ the 
Sarma state exists almost in the region a; > 1 where the 
BCS solution y = \ disappears. Obviously, in a narrow 
region x < 1 there exists a branch of the conventional 
type which is unstable, and the multi-value behavior of 
y means a first order BCS-Sarma phase transition at a 
critical field a;i which is slightly smaller than 1. 

To discuss the thermodynamic stability of the Sarma 
state, we then need to compare it with the normal state. 
In the case of A2 < /i, we find 



No 



2h^ - A^ - 2h\/h^ - Al 



(35) 



Some analytical estimations can be made. For large 
asymmetry A2 <C Ai, around the /i- window h ^ A20 
but h <C Aio for the Sarma state, the sign of the quan- 
tity f2s — r^N is dominated by the first term, if A^i is not 



< 0.6 - 
< 



/(C) 




x=h/A„ 



20 



FIG. 4: The solution y = A2/A20 to the gap equations as 
a function of x = /i/A2o. The dashed line denotes the BCS 
solution y = 1 in the region < x < 1. The solid lines (a) and 
(b) are the Sarma solutions for J 7^ 0. In the calculations, 
we take NxjN-i = 1.5 and JI^/ThUl = 0.07. For (a) we take 
{lJiNx)l{\J2'N2) = 3 which leads to A10/A20 =^ 9, and for 
(b) we have (f/iiVi)/(l72A'"2) = 1.67 and hence A10/A20 ^ 3. 
The dot-dashed line (c) is the conventional Sarma solution for 
J : 



0, 



y 



\J2x — 1 in the region 0.5 < x <\. 



much smaller than A''2. In this case, the BCS solution is 
absent and the Sarma state is the stable ground state. 
This argument confirms our conclusion from the numer- 
ical results in Figl2] and FigO The h window for stable 
Sarma state is wider when the asymmetry between the 
two bands becomes larger. This means that the CC limit 
of such a two-band superconductor can be much higher 
than the conventional value h^. = 0.707A20. 



IV. SUMMARY 

We have studied the stability of Sarma state in two- 
band Fermi systems via both the stability analysis and 
analytical solution at weak coupling. From the stabil- 
ity analysis, the Sarma state can be the minimum of the 
thermodynamic potential and hence a possible candidate 
of the ground state, if the inter-band exchange interac- 
tion is turned on. Both numerical and analytical studies 
show that, a large asymmetry between the two bands 
or the two pairing gaps is an important condition for 
thermodynamic stability of the Sarma state. When the 
condition is satisfied, two first order phase transitions 
will occur when the mismatch increases, one is from the 
BCS to Sarma state at a lower mismatch and the other 
is from the Sarma to normal state at a higher mismatch. 
Our predictions could be tested in multi-band supercon- 
ductors and ultracold atom gases, and such a gapless 
superconductor may have many unusual properties, such 
as magnetism and large spin susceptibility^li^^. 

Acknowledgments: We thank W.V.Liu and M.Iskin 
for useful communications and Y.Liu and X.Hao for the 



8 



help in numerical calculations. The work is supported National Research Program Grant 2006CB921404. 
by the NSFC Grants 10575058 and 10735040 and the 



1 G.Sarma, J.Phys.Chem.SoIid 24,1029(1963) 

2 P.Fuldc and R.A.Ferrell, Phys. Rev. A135, 550(1964); 
A.I.Larkin and Yu.N.Ovchinnikov, Sov.Phys. JETP 20, 
762(1965) 

^ B.S.Chandrasekhar, Appl. Phys. Lett. 1,7(1962); 

A. M.Clogston, Phys. Rev. Lett. 9, 266(1962) 

* K. Gloos, R.Modler, H.Schimanski, C.D.Brcdl, C.Gcibcl, 
F.Stcglich, A.I.Buzdin, N.Sato and T.Koniatsubara, 
Phys. Rev. Lett. 70, 501 (1993); R. Modlor, 
P.Gegenwart, M.Lang, M.Deppe, M.Weiden, T. Lh- 
mann, C.Geibel, F.Steglich, C.Paulsen, J.L.Tholence, 
N.Sato, T.Komatsubara, Y.nuki, M.Tachiki and 
S.Takahashi, Phys. Rev. Lett. 76, 1292(1996); A.Bianchi, 
R.Movshovich, C.Capan, P. G.Pagliuso and J.L.Sarrao, 
Phys. Rev. Lett. 91 (2003) 187004 

^ J. L. O'Brien, H. Nakagawa, A. S. Dzurak, R. G. Clark, 

B. E. Kane, N. E. Lumpkin, R. P. Starrett, N. Muira, E. 
E. Mitcheh, J. D. Goettee, D. G. Rickel, and J. S. Brooks, 
Phys. Rev. B 61, 1584 (2000) 

L. Balicas, J.S.Brooks, K.Storr, S.Uji, M.Tokumoto, 
H.Tanaka, H. Kobayashi, A.Kobayashi, V.Barzykin and 
L.P.Gor'kov, Phys. Rev. Lett. 87, 067002(2001); M. A. 
Tanatar, T. Ishiguro, H. Tanaka and H. Kobayashi, Phys. 
Rev. B 66, 134503 (2002) 
^ M.W.Zwicrlcin, A.Schirotzek, C.H. Schunck and W. Ket- 
terle. Science 311, 492(2006); G.B. Partridge, W.Li, 
R.I.Kamar, Y.-an Liao, R.G. Hulet , ibid.311, 503(2006) 

* P.F.Bedaque, H. Caldas and G. Rupak, Phys. Rev. Lett. 
91, 247002(2003); A.Sedrakian, J. Mur-Petit, A. Polls 
and H. Miither, Phys. Rev. A72, 013613(2005); L.He, 
M.Jin and P.Zlmang, Phys. Rev. B73, 214527(2006); 
Phys.Rev.B74, 214516(2006) 

^ C.H.Pao, S.-T. Wu and S.-K.Yip, Phys. Rev. B73, 
132506(2006); D.E.Sheehy and L.Radzihovsky, Phys. Rev. 
Lett. 96, 060401(2006); D.T.Son and M.A.Stephanov, 
Phys. Rev. A74, 013614(2006); H.Hu and X.Liu, ibid.73, 
051603(R)(2006); M.Mannarelh, G. Nardulh and M. Rug- 
gicri, tbid.74, 033606(2006); A.Bulgac and M.M.Forbes, 
ibid.75, 031605(R)(2007); 

" J.Kinnunen, L.M.Jensen and P.Torma, Phys. Rev. 
Lett. 96, 110403(2006); P.Pieri and G.C.Strinati, 
ibid.96, 150404(2006); T.Silva and E.Mueller, md.97, 
070402(2006); T.Mizhushima, T.Mizushima and M. 
Ichioka,i6id97, 120407(2006) 

^ D.E.Sheehy and L.Radzihovsky, Ann. Phys. (N.Y.)322, 



1790(2007) 

A.J.Leggett, in Modern trends in the theory of condensed 

matter, Springer- Verlag, Berlin, 1980, pp. 13-27 
" M.Huang, P.Zhuang, and W.Chao, Phys. Rev. D67, 

065015(2003); I.Shovkovy and M.Huang, Phys. Lett.B564, 

205(2003); M.Alford, C.Kouvaris and K. Rajagopal, Phys. 

Rev. Lett. 92, 222001(2004) 
" D.Ebcrt and K.G.Klimonko, J. Phys.G32, 599(2006); L.He, 

M.Jin and P.Zhuang, Phys.Rev.D74, 036005(2006). 

W.V.Liu and F.Wilczek, Phys. Rev. Lett.90, 047002(2003) 
1® M.M.Forbes, E.Gubankova, W.V.Liu and F.Wilczek, 

Phys. Rev. Lett. 94, 017001(2005) 
" E.Taylor, A. Griffin and Y.Ohashi, Phys. Rev. A76, 023614 

(2007) 

H. Suhl, B.T.Matthias and L.R.Walker, Phys. Rev. Lett.3, 
552(1959) 

V.M.Loktev, R. M. Quick and S. G. Sharapov, Phys. Rept. 
349,1(2001) 

^° M. lavaronc, G. Karapotrov, A. E. Koshelev, W. K. Kwok, 

G. W. Crabtrcc, D. G. Hinks, W. N. Rang, E.-M Choi, 

H. J. Kim, H.-J. Kim, and S. I. Lee, Phys. Rev. Lett. 89, 
187002 (2001); F. Bouquet, Y. Wang, I. Sheikin, T. Plack- 
owski, A. Junod, S. Lee and S. Tajima, ibid.89, 257001 
(2002); S.Tsuda, T. Yokoya, Y. Takano, H. Kito, A. Mat- 
sushita, F. Yin, J. Itoh, H. Harima and S. Shin, Md.91, 
127001 (2003); J. Geerk, R. Schneider, G. Linker, A. G. 
Zaitsev, R. Held, K.-P. Bohnen and H. v. Loehneysen, 
ibid.94, 227005 (2005) 

^1 M. Kohl, H. Moritz, T. Stoferle, K. Giinter, and T. 

Esslingcr, Phys. Rev. Lett. 94, 080403(2005); L.M.Duan, 

Phys. Rev. Lett.95, 243202 (2005) 
^2 M. Iskin and C. A. R. Sa de Melo, Phys. Rev.B72, 024512 

(2005); ibid.74, 144517(2006) 

J.-P. Martikainen, E. Lundh and T. Paananen, Phys. 
Rev.A78, 023607(2008) 

A.I.Akhiezer, A.A.Isayov, S.V.Polctminsky and A. A. Yat- 
senko, Phys. Lett. B451, 430(1999) 

R.Casalbuoni and G.NarduUi, Rev. Mod. Phys. 76, 

263(2004) 

^'^ D.E.Sheehy and L.Radzihovsky, Phys. Rev. B75, 
136501(2007) 

S.Takada and T.Izuyama, Prog.Theor.Phys.41,635(1969). 
2* L.He, M.Jin and P.Zhuang, Phys.Rev.B73, 024511(2006) 



